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AN ASYMPTOTIC SHAPE THEOREM FOR RANDOM LINEAR 

GROWTH MODELS 

AURELIA DESHAYES 


Abstract. In this note, we generalize the asymptotic shape theorem proved in |DesI4a| 
for a class of random growth models whose growth is at least and at most linear. In 
this way, we obtain asymptotic shape theorems conjectured for several models: the 
contact process in a randomly evolving environment |SW08| . the oriented percolation 
with hostile immigration |GM12b| and the bounded modified contact process |DS00| . 


1. Introduction 

In 1974, Harris |Har74| introduced the classical contact process as an interacting parti¬ 
cle system modeling the spread of a population through the grid This process has been 
well studied and it satisfies in particular a shape theorem. Since then, a lot of extensions 
appeared in the literature: two stage eontact process by Krone |Kro99) . boundary modi¬ 
fied contact process by Durett and Schinazi |DS00| . contact process in randomly evolving 
environment by Broman |Bro07| . Remenik |R,em08) . Steif and Warfheimer jSWOS], con¬ 
tact process with aging by the author |Desl4a| etc. All these processes are linear random 
growth models in the sense that we can exhibit a quantity of interest At taking values in 
V{7f) and deterministic compact sets A_ and A_|_ such that, for t large enough, 

C Ai C tA+. 

The first inclusion traduces the at least linear growth and the second one the at most 
linear growth. These models have also others characteristics like attractivity, phase tran¬ 
sition phenomenon etc. Here, we want to show that they satisfy asymptotic shape results 
i.e. there exists a deterministic set B such that ^ converges to R in a sense to precise 
and we want to highlight their common properties which lead to a general asymptotic 
shape theorem. In order to do that we introduce a general class of random growth models 
including the above mentioned contact process extensions and a related good quantity, 
the essential hitting time a, which allows us to use the same techniques that the ones 
involved in |Desl4a| to prove an asymptotic shape theorem. Indeed, as the standard con¬ 
tact process itself, its extensions are non permanent models: the extinction is possible 
and the hitting times can be infinite so that the standard integrability conditions are 
not satisfied. On the other hand, if we condition the model to survive, then stationarity 
and subadditivity properties can be lost. The quantity a overcome such lacks: it turns 
out that this function a satisfies adequate stationarity properties as well as the almost- 
subaddivity conditions involved in Kesten and Hammersley’s theorem (see |Kes73) and 
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|Ham74) '). a well-known extension of Kingman’s seminal result. We will define precisely 
all the quantities involved but we will refer to |Desl4a) for proofs working in the same 
way. 

We obtain a general asymptotic shape theorem for a class of random linear growth 
models and deduce asymptotic shape theorems conjectured for several models: the con¬ 
tact process in a randomly evolving environment |SW08) . the oriented percolation with 
hostile immigration [GM12b) and the bounded modihed contact process [DSOD) . 


2. A CLASS OF RANDOM LINEAR GROWTH MODELS 

2.1. Notations and definitions. 


2.1.1. State space. We will work on interacting particle systems on For x and y in 
Z<^, we say that x and y are neighbors, and we denote by x ~ y, if 1^* ~ Vil ~ 

We denote by Tx the spatial translation operator of vector x £ Z'^. 

Let S' be a finite totally ordered set which will represent the possible state (or type) 
of a particle x £ Z'^. We denote by S^ the set of mappings ^ : Z”* i-A S which can also 
be seen as the set of partitions of Z'^ in \S\ subsets: 

Z^ ^ S 






X 


and we denote by Sj the subset of functions with finitely many coordinates non equal 
to minS (often 0). Let be the set of continuous functions from S^'^ to M and 

Co(S^'*) the subset of functions depending on only hnitely many coordinates of The 
set will be our state space of configurations. 


2.1.2. Markov process and property of interest. We work with a time space T which can 
be discrete (T = N) or continuous (T = K+)- Let {^t)teT be a stationary Markov process 
taking values in and let (P) be a property of interest about elements of S. We define 
the process (At) taking values in V{Z'^) by, for t £T, 

At = {x £ Z'^,f^t{x) satisfies (P)}; 

the quantity At represents the set of points which satisfies the property (P) at time t. 

We can also write, for t £ T and x £ Z'^, At{x) = g{f,t{x)) where g is the indicator 
function corresponding to the property (P); we extend At = g{ft) in the natural way. 
One should notice that the process {At) is not necessarily Markovian. 

Let (imin be the minimal non trivial configuration of the process {ft) such that Aq = 
g{Smm) is not empty {i.e. someone in (imin satisfies (P)). A such configuration (5min exists 
because S is totally ordered and it is unique modulo translation. For every x £ Z^, we 
denote by 

r = inf{t > 0, = 0}, 

t{x) = inf{t > 0,x G Aj™}. 

Starting from (iminj x is the first time when no one satisfies (P) and t{x) is the first time 
when X satisfies (P). 
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For X G and t gT we denote by Tx the spatial translation and by 9t the temporal 
one. For t gT, let 

Ht = [_} At = {x <t} 

S<t 

be the set of points which have satisfied (P) before time t. 

2.1.3. Interacting particle system. 

Definition. Let {^t) be a Markov process taking values in and (P) a property of 
interest. We say that belongs to the class C if it satisfies the following conditions: 

(1) The generator A of (f^t) satisfies, for every f G Co{S^'^) and every configuration 
f G 

Af)iO = '^c{x,C,s) {f{Cix,s)) - fiO) 

x&'i ses 

where, for all x G lA and ^{x,s) tb.e configuration ^ where ^(x) is replaced 
by s. c(x, s) >0 is the intensity for a jump ^ depends on ^ only 

through {£^{y),y ^ x} and ^{x) and we take c{xA,s) = Q if s = f{x). 

(2) The process is additive in the following sense: for all f,g G = 

iti ‘moreover, {At = 0 } in an absorbing state. 

The first condition implies that is a nearest neighbor interacting particle system; 
this definition was introduced by Harris in |Har74| but today the term NNl represents a 
larger class of interacting particle systems (including for example the exclusion process 
which does not satisfy this first condition). We can find the general definition of NNl 
particle systems and the proof of their existence in |Lig05| . 

Harris [Har78) introduced a graphical representation of additive processes which make 
possible to work with a percolation view on a diagram of Poisson processes impacts 
(see |Desl4a| for an example of such construction). 

2.1.4. Random linear growth models. We define the class of models with which we are 
going to work: 

Definition. We say that belongs to the class Cl if it satisfies the following con¬ 

ditions: 

(1) we have P(r = +oo) > 0 

(2) there exist Ci, C 2 , Mi, M 2 > 0 such that for all t > 0 and x GTA 

(AML) P ^3x € : t{x) < t and ||x||i > Mit^ < Ci exp(— 6 * 2 ^), 

(SC) P (t < r < 00 ) < Cl exp(— 6 * 2 ^), 

(ALL) P (t(x) > M 2 ||x|| + t, r = 00 ) < Ci exp{—C 2 t). 

The second condition means that the growth of the set of points satisfying (P) is 
at most linear (AML), at least linear (ALL) and, if the extinction time r is finite. 
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then it is small (SC for small clusters in analogy with percolation vocabulary). Since 
P(r = +oo) > 0, we will work with the probability conditioned to survive: P = P(.|t = +oo). 

2.2. Examples. 

2.2.1. Classical contact process. For x G ^ G {0,1}^"^ and /c G N: 

C(x,^,0) = l^(a;)=l 

We consider P the property ’to be alive’ and we get back to the already known asymptotic 
shape theorem for At = it- In the following we denote by Ac the critical value for classical 
contact process. 

Next examples are extensions of the classical contact process from the literature; we 
give here dehnitions with the notations of our context. They belong to class C by con¬ 
struction. Several authors have worked hard to show that they belong to the class Cl- 

2.2.2. Contact process in a randomly evolving environment (CPREE). The contact pro¬ 
cess in a randomly evolving environment was introduced by Broman |Bro07) and studied 

in particular by Steif and Warfheimer |SW08) . Take T = M+, S = {(0,0), (0,1), (1,0), (1,1)}. 
An element of S represents a pair (type, stage) where type can be favorable ( 1 ) or unfa¬ 
vorable (0) and the stage can be alive (1) or dead (0). We can write the jump intensities: 



i[x) ^ s 

c{x,i,s) 

birth 

( 0 , 0 ) ^( 0 , 1 ) 


(type independent) 

( 1 , 0 ) ^( 1 , 1 ) 

'^x^v ^t{y) 

death 

( 0 , 1 ) ^( 0 , 0 ) 

So 

(type dependent) 

( 1 , 1 ) ^( 1 , 0 ) 

Si 


( 0 , 0 ) ^( 1 , 0 ) 

IP 

evolution 

( 0 , 1 ) ^( 1 , 1 ) 

IP 

of type 

( 1 , 0 ) ^( 0 , 0 ) 

7(1 -p) 


( 1 , 1 ) ^( 0 , 1 ) 

7(1 -p) 


with 0 < < 5q, p ^ [ 0 , 1 ] and 7 > 0 four parameters of the system. 

Let (P) =’to be alive’. The process [At) represents the set of alive points regardless 
of their type (favorable or unfavorable) and corresponds to the process [Ct) in Steif and 
Warflieimer’s article [SW08) . Here, (imin is the conhguration where 0 is in state (0,1) and 
all the others are in state (0,0). Using the Steif-Warfheimer construction and a restart 
argument like the one in section 5.3 in [Desl4a| . it is easy to see that the controls (AML), 
(SC) and (ALL) are satished, so [it,P) belongs to C CCl- 

2.2.3. Contact process with aging (CPA). We introduced the contact process with aging 
in [Desl4a) where the dynamics are detailed; in this article we made a precise construction 
which allows us to consider a general case where a particle has an integer age. In the 
present context we have to suppose S finite so we can consider an aging contact process 
with a maximal age N. Then take T = M_|_ and S = {0,..., A^}. 


















AN ASYMPTOTIC SHAPE THEOREM FOR RANDOM LINEAR GROWTH MODELS 


5 


Let (P) =’to be alive’. The process [At) represents the set of alive points regardless of 
their age. Here, Jmin = <5o is the configuration where all particles are dead except 0 which 
is alive with age 1. From the results obtained in [Desl4a| . belongs to C nC^. 

2.2.4. Dependent oriented percolation (DOP). In |GM12b| . Garet and Marchand intro¬ 
duced a dependent oriented percolation model where fertile bacterium (represented by 
a type 1 particle) is submerged in a population of immune cells (type 2 particle) that 
are going to impede its development. The immune cells are not very fertile but benefit 
from a constant immigration process. Take T = N and S = {0,1,2}. The system is 
described by a discrete time Markov chain depending on 3 parameters p,q,a € (0,1). 
Firstly, between time n and time n+ 1/2, each particle attempts to infect its neighbors: 
it succeeds with probability p if it is a type 1 particle, and with probability q if it is 
a type 2 particle (squeezing a type 1 particle if there is one). Secondly, between time 
n -|- 1/2 and time re -|- 1, there is an immigration of type 2 particles on each site with 
probability a > 0 (again eventually squeezing a type 1 particle). Let (P) =’to be of 
type 1’. The process {At), denote by (iyi,t) in |GM12bj . represents alive fertile particles. 
Here, dmin = where all particles are dead except 0 which is alive and fertile (of type 1). 
Theorem 1.2 of [GM12b) assure that the controls (AML), (SC) and (ALL) are satisfied if 
the process survives (that is if a is smaller than a non trivial value ac{p,q))- So, (^t,P) 
belongs to C G 


2.2.5. Boundary modified contact process (BMCP). Durrett and Schinazi [DS00| intro¬ 
duced a Boundary Modified Contact Process where particles can be infected, susceptible 
that has never been infected (state —1) and susceptible that as been previously infected. 
Take T = M+, S = { — 1,0,1} and the jump intensities as follows: 


c{x,r], 1) = 


and c{x,rj,0) = 1 if r/(x) = 1. 


\eNi{x,fi) if 77(x) =-1, 

^\iNi{x,Tj) if ry(x) = 0, 

The quantities Ag and Aj are two non negative parameters of the system. The process 
[At] represents infected particles at time t. Here, dmin is the la configuration where all site 
are susceptible never infected {i.e. in state —1) except site 0 which is infected (state 1). 
We think that the renormalization work done by Durrett and Schinazi in |DS00| makes 
possible to obtain the controls (AML), (SC and (ALL) thanks to restart techniques (like 
in section 5.3 of |Desl4a| l for A* > max(Ac, Ag). 

For contact processes (with aging, boundary modified, in a randomly evolving envi¬ 
ronment), t{x) is the first time when the particle x is alive regardless its age, type or 
memory. For dependent oriented percolation, t{x) is the first time when the particle x is 
one of type 1 alive 


2.3. Results. In this note, we give the elements to show the following general asymptotic 
shape theorem. 

Theorem 1. Let {^t) be a Markov process taking values in and (P) an associated 
property of interest. If {^t,P) belongs to classes C and Cl then there exists a norm pL on 
such that for all e > 0, P(.|r = -|-oo) almost surely, for all t large enough, 

(1 — e)B^ C — C (1 -I- e)R^, 
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with Gt = {x ^ Z'^,t(x) < t} + [0, l]'^ and the unit ball associated to fi. 

We can immediately deduce asymptotic shape theorems for the above examples. In 
each statement, we denote by the unit ball associated to the norm /i and P the 
probability conditioned to the survival of the property (P) starting from the configuration 

First we recall theorems already proved. 

Theorem 2 (Asymptotic shape theorem for CPA). Let (^t) a contact process with aging. 
If {^t) is supercritical, then there exists a norm pL on such that for all e > 0, P almost 
surely, and for all t large enough 

(1 — £)B^ C — C (1 + £)B^, 

where Ht = Ht + [0,1]'^ with Ht the set of points born before t, regardless of their age. 

We deduce from Theorem [2] an asymptotic shape theorem for Krone’s model [Kro99| 
which is a particular case of the contact process with aging. 

Now we announce theorems previously conjectured. 

Theorem 3 (Asymptotic shape theorem for CPREE). Let (ft) a contact process in a 
randomly evolving environment. If p > pc (critical value for bond percolation on 
then there exists a norm p on such that for a// e > 0, P ( . = +oo) almost 

surely, for all t large enough, 

(1 — £)B^ C — C (1 + £)Bfj_, 

where Ht = Ht + [0,1]'’* with Ht the set of points born before t, regardless of their type. 

Theorem 4 (Asymptotic shape theorem for DOP). Let (rjt) a dependent oriented per¬ 
colation. If p > pi, (critical value for oriented bond percolation on x 'N), q < pi 
and a G {0,ac{p,q)), then there exists a norm p on such that for all £ > 0, 
P( . = +oo) almost surely, for all t large enough, 

(1 — £)B^ C — C (1 + £)Bfj_, 

where Ht = Ht [0,1]'’* with Ht the set of sites where particles of type 1 are before t. 

Theorem 5 (Asymptotic shape theorem for BMCP). Let (pt) a contact process boundary 
modified. If A* > max(Ac, Ag) and subject to that there exist Ci,C 2 , Mi, M 2 > 0 such that 
for all t > 0 and x £ 

P G Z^ : t{x) < t and ||x||i > Mit^ < Ci exp(—C' 2 t), 

P (t < r < 00 ) < Cl exp(—C' 2 t), 

P {t{x) > M 2 ||x|| +t,T = 00 ) < Cl exp(—C' 2 t), 

oil T = inf{t > 0,Vx G Z'^,pt{x) < 0} et t{x) = inf{t > 0,r/t(x) = 1}, then there exists 
a norm p on such that for all e > 0, P( . |r = + 00 ) almost surely, for all t large 
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enough, 

(1 — e)B^ C — C (1 + 

where Ht = Ht + [0,1]“^ with Ht the set of particles of type 1 born before t. 

From now, (^t, P) belongs to CCiCl- In Section[H]we introduce the essential hitting time, 
a quantity that has good properties regarding to the dynamical system. In Section U] we 
prove the Theorem [1] thanks subadditive ergodic theorems applied to the essential hitting 
time. 


3. Essential hitting time 


3.1. Definition. With non permanent models like contact process extensions, the hitting 
times can be infinite (because extinction is possible) and if we condition on the survival, 
we can lose independence, stationarity and even subadditivity properties required by 
Kingman theory. 

We are inspired by the construction of Caret and Marchand |GM12a| for contact 
process in random environment to build the essential hitting time for a model which 
belongs to C PI Cl . 

We set uq{x) = Vo{x) = 0 and we define by induction two sequences of stopping times 
{un{x))n and {vnix))n as follows. 

• Suppose that Vk{x) is defined. We set 

Uk+iix) = inf{t >Vk{x) : X G 

If Vk{x) is finite, then Uk+i{x) is the first time after Vk{x) where the site x satisfies 
(P); otherwise, Uk+i{x) = +oo. 

• Suppose that Uk{x) is defined, with A: > 1. We set Vk{x) = Uk{x) + o 9u^,{x)- 
If Uk{x) is finite, then the time r* o 9uf,{x) is the (possibly infinite) extinction 
time of (P) starting at time Uk{x) from the configuration (fmin o Tx] otherwise, 
Vk{x) = + 00 . 

We have uo{x) = vo{x) < ui{x) < vi{x) < ... < Ui{x) < Vi{x) ... We then define K{x) 
to be the first step when Vk or Uk+i becomes infinite: 

K{x) = min{A: > 0 : Vk{x) = +oo or Uk+i{x) = +oo}. 

K{x) has a sub-geometric tail. In particular, K(x) is almost surely finite. It allows us 
to define the following quantities. 


Definition. For x we call essential hitting time of x by the property (P) the 

quantity cr{x) = Uj^^x)- define the operator 6x on by setting: 


( 1 ) 


Ox 


Txo0^(^x) if(r{x)<+oo, 
Tx otherwise. 


Thanks to controls (AML), (SC) and (ALL), we can obtain that, for x G \ {0}, P 
is invariant under the action of Ox- Besides, under P, <j{y) o Ox is independent from <t(x) 
and its law is the same as the law of CF{y). We refer to [Desl4a) for proofs of the previous 
properties about K{x) and Ox (and to chapter 7 of |Desl4b| for more details). 
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3.2. Bad growth points. In order to apply the techniques used in Section 7 of |Desl4a| 
to control a, we introduce the general box of bad growth point for a model in C H 
For every x G L > 0 and t > 0, we introduce the event 

Ey{x, t) = {Hf <^y + BMxt} U {t/2 < < +oo} U {r^ = +oo, inf{s > 2t : x G A^} > Kt} 

traducing the fact that a point y has bad growth with respect to the spatio-temporal point 
(x,t). Let Vy the counting measure of all events occurring at y during [0, L]. Thanks to 
Harris representation, can be expressed with Poisson processes. This measure allows 
us to count the number of bad growth points by 

NL{x,t)= 

y£x+BMt+2 

For technical reasons, we add 0 in the measure u and we modify Ey{x,t) by 

( 2 ) Ey{x,t) = Ey{x,t)^{i'y[0,t/^] = 0 } 

Going back to our examples: 

• For DOP (which is a discrete model), u is the counting measure on N. 

• For CPA we consider: 

y'y — E + 5o 

eGM^-.yGe 

where iVy, cOy and are the Poisson processes respectively giving the possible 
death times at y, the possible maturations and the possible birth times through 

e. 

• For CPREE, we consider: 

^y — ‘^y E UJy + Ue + 5o) 

eGE'^:y€e 

where ojy is the Poisson process associated to type 0 deaths, Uy the one associated 
to type 1 deaths and ojf. the one associated to births. 

With this definitions we can show the following properties. 

Lemma 6. Let x G Z'^ and t > 2. If L,s are positive integers such that NL{x,t) o9s = 0 
and Ui{x) G [s -|- t, s -|- L], then Vi{x) = -|-oo ou ttj+i(x) — Ui{x) < nt. 

Besides, for any t > 0, s G N and x G the following inclusion holds: 

{t = -l-oo} 

(3) n {3u G X + BMit+ 2 -, o 0s = -|-oo, u G As} 

(4) n {NK(x)Kt{x, t) o 0s = 0} 

(5) n Pi {vi{x) - Ui{x) < t} 

l<i<K (fc) 

C {r = -|-oo} n {(t(x) < s -t- K{x)Kt} . 

Then, following the same proof lines that |Desl4a) . and choosing k = 3Mi(l -|- M 2 ), 
where Mi and M 2 are the constants respectively given in (AML) and (ALL), we control 
the probability that a space-time box contains no bad growth point. 
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Lemma 7. There exist A,B>{) such that for all L > 0, x £ and t > 0 one has 

P(A/i(x,i) > 1) < Aexp{—Bt). 

3.3. Subadditivity and difference between a and t. In the same way, the following 
controls are still available. 

Proposition 8. (I) There exists A,B>0 such that for all x,y £ 

Vt > 0 P + y) “ + ^{y) ° < ^exp (^—By/t^ . 

Moreover, for p> 1, there exists Mp > 0 such that for all x,y £ Z'^, 

E[((t(x + y) - (cr(x) + a{y) o < Mp. 

_ — i(x}| 

(II) almost surely, it holds that lim -r — r- = 0. 

||3l||^ + 00 ||x|| 

(III) There exist A,B,C > 0 such that 

Vx £ Z'^,Vt > 0, P((t(x) > Cllxll +t) < Aexp{B^/t). 

Moreover, for p > 1, there exists a constant C{p) >0 such that 
Vx E Z^ E[(j(xf] < C{p){l + ||x||f. 

(IV) For every e > 0, P-a.s., there exists R > 0 such that 

Vx,y £ Z'^, (||x|| > R and ||x — y\\ < e||x||) (|o'(x) — cr{y)\ < (7611x11) . 


4. Proof of the asymptotic shape theorem 


Thanks to the almost subadditive theorem of Kesten and Hammersley (see |Ham74) 
and |Kes 73]), we obtained a general asymptotic shape theorem on some random variables 


Theorem 9 (Theorem 39 of |Desl4a) l. Let (fl, R, P) be a probability space. Let {(y{x))^^i^d 
be random variables with finite second moments and such that, for every x £ , o'(x) 

and o'(—x) have the same distribution. Let {s{y))y^id and {r{x,y))^ y^j^d be collections 
of random variables such that: 

Hyp 1: Vx,y E Z*^, cr(x + y) < <t(x) + s{y) + r{x,y) with s{y) having the same law as 
cr{y), and being independent from cr{x), 

Hyp 2: Vx, y £ Z'^, ^Cx,y and ax,y < 2 such that 
Vn,p,E[r{nx,py)‘^] < Cx,y{n + p)'^^'y, 

Hyp 3: 3(7 > 0 such that Vx E Z*^, P(cr(nx) > (7n||x||) 0, 

Hyp 4-' 3iV > 0 such that Ve > 0,P — p.s 3M such that (||x|| > M and \\x — y\\ < 
K\\x\\) => ||o'(a;) - (j{y)\\ < e||x||, 

Hyp 5: 3c > 0 such that Vx E Z*^, P(cr(nx) < cn||x||) 0. 

Then there exists /r : Z'^ —>■ M'*' such that 


lim 

||a:i||^oo 


O'(x) — p{x) 


= 0 a.s. 


X 












10 


AURELIA DESHAYES 


Moreover, fi can be extended to a norm on and we have the following asymptotic shape 
theorem: for all e > 0 , P almost surely, for all t large enough, 

(1 — e)B^ C — C (1 + e)B^, 

where Gt = {x G : (t{x) < t} + [ 0 , 1 ]'^ and B^ is the unit hall for p,. 

We now deduce the expected asymptotic shape theorem for the hitting time t: 


Proposition 10. 


There exists a norm p, on R*^ sueh that almost surely under P, 

lim ^0, 

||fc||—>-+oo ||x|| 


and for every e > 0, P-a.s., for every large t, 

(1 - e)Bfj, C C {1 + €)Bfj_ 

where = {x G Z'^/t(x) < t} + [0, l]'^ and B^ is the unit ball for p. 
Theorem [1] is contained is the previous result. 


Proof. First, we use Theorem[9]to show that a satisfies an asymptotic shape theorem. We 
check the hypotheses of Theorem [9] using the controls of Proposition | 8 j Thanks to m, 
a has finite second moment required. We take s(y) = cr(i/) o 0^. The hypotheses 1 and 
2 are satisfied thanks to properties of Ox and ([T|). The hypothesis 3 is the at least linear 
growth m and the hypothesis 5 is immediately checked thanks to the at most linear 
growth (AML): P(iT(nx) < M 2 n||x||) < P(t(nx) < M 2 n||x||) < ^ exp(—SM 2 n||x||). 
Finally, the hypothesis 4 is the control m- We deduce the result for t from the result 
for a thanks to (P). □ 
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